Quantum-dot Carnot engine at maximum power 
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We evaluate the efficiency at maximum power of a quantum-dot Carnot heat engine. The universal 
value of the coefficients at the linear and quadratic order in the temperature gradient are reproduced. 
Curzon-Ahlborn efficiency is recovered in the limit of weak dissipation. 



I. INTRODUCTION 



The purpose of a heat engine is to transform an amount 
of heat Qh, extracted from a hot reservoir at temper- 
ature T/j, into an amount of work W. The efficiency 
j] = W/Qh for doing so is at most equal to Carnot ef- 
ficiency: r\ < r) c , with r\ c — 1 — T c /Th- Here T c is the 
temperature of a second, cold reservoir T c < X^, in which 
the remaining energy Qh — W is deposited. The equality 
is reached for reversible operation, implying that the cor- 
responding power output is zero. Curzon and Ahlborn 
were amongst the first to study the question of efficiency 
at maximum power By considering a simple modi- 
fication of the Carnot engine and after applying the so- 
called endo-reversible approximation (neglecting dissipa- 
tion in the auxiliary system), they found the following 
beautiful expression for the efficiency at maximum power: 
Vca = 1 - VTjn = 1 - v^T - Ifc '= r lc /2 + 7 ?c 2 /8 + • • • . 
While this formula appears to describe rather well the 
efficiency of actual thermal plants, and is close to the 
efficiency at maximum power for several model systems, 
it is neither an exact nor a universal result, and it is 
neither an upper nor a lower bound. It has therefore 
come as a surprise that a number of universal predic- 
tions can be made about the expansion of the efficiency 
at maximum power in terms of rj c . In the regime of lin- 
ear response, i.e., at first order in r\ c , it is found that 
the efficiency at maximum power is at most half of the 
Carnot efficiency 0). In other words, the CA efficiency is 
an upper bound at the level of linear response. The proof 
was given for systems operating under steady-state condi- 
tions. The upper bound is reached for so-called strongly 
coupled systems, i.e., systems in which the heat flux and 
the work-producing flux are proportional to each other. 
More recently, it has been shown that the quadratic co- 
efficient, equal to 1/8, is also universal for strongly cou- 



pled systems in the presence of an additional left-right 
symmetry in the system [3J. Furthermore, the univer- 
sality of the coefficients is a direct consequence of the 
time-reversibility of the underlying physical laws. The 
coefficient 1/2 derives from the symmetry of the Onsager 
matrix. The coefficient 1 /8 can be seen as the implication 
of Onsager symmetry at the level of nonlinear response. 

The above universality predictions have been con- 
firmed in a number of steady-state model systems in- 
volving classical particles [|[ , fermions @ , and bosons [f| . 
Universality has also been observed in variants of the CA 
model based on Carnot cycles performed in finite time, 
even though finding the optimal driving protocol max- 
imizing work extraction can be notoriously difficult 0- 
[l3j . The connection with the steady-state analysis has 
been clarified by identifying the Onsager coefficient for 
a finite-time Carnot cycle in the linear regime 11]. Fur- 
thermore, agreement with the universal quadratic coeffi- 
cient has also been observed in a Carnot cycle based on 
a (classical) Brownian particle in an harmonic trap @. 

In this paper we provide a full analysis of a thermal en- 
gine undergoing a Carnot cycle, with the auxiliary system 
consisting of a single-level quantum dot that is switched 
between a hot and a cold reservoir. We show that the ef- 
ficiency at maximum power is again consistent with the 
above discussed universality. Furthermore, CA efficiency 
at maximum power is obtained exactly at all orders in 
■q c in the limit of low dissipation, which is distinct but 
similar to the case considered in the original CA paper. 



II. MODEL 

Our heat engine model consists of a single-level quan- 
tum dot interacting with a metallic lead through a tun- 
neling junction. The quantum dot is assumed to have a 
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single energy level e near the Fermi level of the lead while 
other levels in the dot do not contribute to the processes 
described below. The state of the system is specified by 
the occupation probability p(t) of having an electron in 
the dot. The lead plays the role of a thermal bath at 
temperature T and chemical potential fi. Electrons are 
assumed to thermalize instantaneously upon tunnelling 
into the lead. 

When the energy level e is modulated by an external 
agent according to a given protocol, a certain amount of 
energy, positive or negative, flows into the system in the 
form of work and/or heat. In the case of an occupied 
level, an amount of work equal to (e/ — [if) — (ej — /Xi) 
is delivered to the system, where the subscripts / and i 
refer to final and initial values. When the electrons at 
energy level e tunnel in (out), an amount of heat equal 
to Q = e — /i (Q = — e + /i) is extracted from the bath. 

The basic problem that we address is the finite-time 
performance of this engine as it runs through the follow- 
ing four standard stages of a Carnot cycle (also see figure 
®: 

I Isothermal process 
The quantum dot is in contact with a cold lead 
at temperature T c and chemical potential /i c . The 
energy level is raised from £q to £\ according to a 
certain protocol during a time interval of duration 
t c . Both work and heat are exchanged during this 
process. 

II Adiabatic process 

The quantum dot is disconnected from the cold 
lead, and the quantum level is shifted from £\ to a 
new level £2- Since the quantum dot is thermody- 
namically isolated, the population of the level does 
not change during this process. Hence, there is no 
heat exchange. However, the change of the energy 
level releases a corresponding amount of work. We 
assume that the operation time of this step is very 
short, in particular negligibly small compared to 
that of the isothermal processes. 

III Isothermal process 

The dot is connected to the hot lead at temperature 
Th and chemical potential fih- The energy level is 
lowered from £2 to £3 based on another protocol 
during a time interval of length . Just as in step 
I, both heat and work are exchanged. 

IV Adiabatic process 

The system is again disconnected from the lead and 
the level is restored from £3 to the initial level £0, 
at the cost of a corresponding amount of work. Af- 
terwards, the dot is reconnected to the cold lead. 
Again, we assume that the operation time of this 
process is negligibly small. 

The above procedure defines one cycle of the thermal 
engine, requiring a total time r c + 77,. The protocols 
in steps I and III must be designed in such a way that 
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FIG. 1: A Carnot cycle of the model heat engine consisting 
of a single-level quantum dot interacting with a lead through 
a tunnleing junction. 



the thermodynamic state of the system, in our case the 
occupation probability p of the quantum level, returns 
to the same initial value after every cycle. Since there is 
no change in occupation probability during the adiabatic 
stages II and IV, the change in occupation probability 
from, say po to p±, during process I, must necessarily 
be compensated by a change back from p\ to po during 
process III. 

The time evolution of the occupation probability p(t) 
for the state of the quantum dot in contact with a lead at 
temperature /3 _1 (k^ = 1) obeys the following quantum 
master equation: 



p{t) = -oj a {t)p(t) + oj b (t)[l - p(t)}, 



(1) 



where the uj a and wj, are transition rates. In the wide- 
band approximation, these rates are given by 



C 



e -/3[ E (t)-M(t)] + 1 

G 

e +/9[e(t)-M(t)] + 1' 



(2a) 
(2b) 



where C is a rate constant. Noting that raising the energy 
level is equivalent to lowering the chemical potential, we 
introduce an effective energy level e = e— fi. The master 
equation (TTJ can now be rewritten as 



P (t) = -Cp(t) + 



c 



1 



(3) 



The effective level varies along the Carnot cycle as eo = 
£0 - Mc ei = £1 - Mc ->• £2 = £2 - Ph £3 = £3 - M/i- 
Note that the change in the chemical potential is included 
in the jump of the effective level during processes II and 
IV. 

We next turn to the thermodynamic description of the 
model. We use the convention that heat entering the 
system is (like work) positive. The internal energy of the 
system at time t is 



where 



£(*) = U(t) -nK{t) = e(t)p(t), 



U(t) = £(t) P (t), 7i{t)=p{t). 



(4) 



(5) 
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The rate of change of the internal energy, £, is the sum 
of two parts, namely, a work flux W and a heat flux Q, 



W = ep = ep — /ip 
Q = ep = ep — [ip. 



(6a) 
(6b) 



Note that the particle exchange contributes to the heat 
flux [last term in Eq. (|6b|)]. When the energy level is 
below the chemical potential, the direction of heat flow 
is opposite to the direction of tunnelling. 

The net total work and net total heat during the pro- 
cess of duration r are obtained as functionals of the oc- 
cupation probability, 



W[p(-)] = / e(t)p(t)dt 
Jo 



e(t)p(t)dt. 



(7a) 
(7b) 



For cyclic processes, we have £(0) = £(t) and hence W + 
Q = 0. For mathematical simplicity, we evaluate power 
using net heat instead of net work: 



-W 

T 



Q 

T 



(8) 



III. OPTIMIZATION: GENERAL CASE 



Our goal is to maximize the power output and to eval- 
uate the corresponding efficiency. Power is a complicated 
functional of the time-dependent protocols in stage I and 
III, and an exact analytical analysis looks difficult at 
first sight. The optimization can however be done in two 
steps. First, we fix parameter values, r c , t%, po and pi 
and maximize the power with respect to the functional 
space of e(t). Since the total operation time is fixed, we 
just need to maximize the heat. Next, we further maxi- 
mize the power with respect to the remaining degrees of 
freedom t c , 77, , p and p\. The problem of maximizing 
heat or minimizing work for a single-level quantum dot 
moving between given initial and final energy states has 
already been analysed in detail in |14| . We reproduce the 
crucial steps of this analysis for self-consistency. 

To find the protocol that maximizes the heat, we do not 
search directly for the optimal schedule e(t), but identify 
the optimal occupation probability p(t). This is done by 
expressing e{t) in terms of p(t) and p(t), and rewriting 
the heat, Eq. (|7b|) . as a functional of p(t) and p{t): 



pQ\p{-)] = / L(p,p)dt, 



where 



L = In 



i 



P(t)+ P (t) 



1 



p(t). 



(9) 



(10) 



The extremum is found via the standard Euler-Lagrange 
method, leading, after integration, to 



P 



dp (Cp + p)[C(l-p)- P ] 



K. 



(11) 



Here K is the constant of integration. Solving the 
quadratic equation for p, we obtain two first order ODEs, 



P = 



K{1 - 2p) =F y/K 2 + AKp(l - p) 
2C(f + K) ' 



(12) 



The upper sign (— ) should be used for upward processes 
in which the quantum level is raised and the lower sign 
(+) for downward processes. It is worth mentioning a 
useful symmetry between electrons and holes. We are 
using the state of an electron, e(t) and p(t), to describe 
the state of system. Instead, we can also use the state 
of holes, — e(t) and 1 — p(t). If p{t) is a solution for an 
upward process, then 1— p(t) is a solution for a downward 
process with — e(t). Hence we do not need to calculate 
the downward process separately, as it follows from this 
symmetry. 

Before turning to the solution of the differential equa- 
tion (|12[) . we examine the physical meaning of the con- 
stant K . Eliminating p in Eq. (jTTJ) by using the master 
equation ([3]), the resulting quadratic equation for p(t) 
leads to the relation 



p(t) 



1 



3 M«) + l 



1 



M*)/2 



(13) 



This equation indicates that when K — 0, p(t) is the equi- 
librium distribution associated with the instantaneous 
value of the energy, implying that K — corresponds 
to the quasi-static limit (r — > 00). As K increases, p(t) 
deviates from the equilibrium distribution. We conclude 
that K measures how far the state of the system deviates 
from the quasi static limit. We will use this insight below 
to obtain a perturbative solution for small dissipation by 
assuming that K is small. 

Next we proceed to solve Eq. (fT2"j) . Separation of the 
variables p and t leads to the following explicit result for 
the upward processes: 



Ct = F[p(t);K}-F[p(0);K}, 



(14) 



where 



F(p-K) 



hip 



1 



arctan 



1 - 2p 



I In 

2 



K [y/K + 4p(l-p) 

2p + K+ y/K 2 +4Kp(l-p) 



2(1 -p) + K + yjK 2 + 4Kp(l - p) 



(15) 



For the downward processes, we need to use F(l—p; K). 

The value of K is determined by the boundary condi- 
tions: 



Ct = F[ P (t);K}-F[ P (0);K}. 



(16) 



Note that K depends solely on the operation time r, the 
probabilities p and p\ , and the tunneling rate C but not 
on temperature. Unfortunately, the function (|15[) is quite 
complicated so we can not obtain an analytical expression 
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for K. In general we need to solve for it numerically. 
However, an exact perturbative solution is possible, cf. 
the next section. 

Having thus obtained the optimal p(t) with K deter- 
mined by (|14[) . we insert this expression in Eq. f|Tb[) to 
obtain the corresponding maximum heat for the optimal 
upward processes, 



which leads to the efficiency of the engine 



/3Q= / e(t)pdt 

p(r) 

dp In 

p(0) 

= §[p(r);K}- 
where 



p(-r) 



p(0) 

2p(l-p) 



s(p)dp 



K + yjK 2 + 4A>(1 - p) 



2p 2 



S[p(0);if] = AS, 



(17) 



S(p; K) = pin 



2(1 -p)p + K - y/K 2 +4Kp(l-p) 



2p 2 



l-2p 



In 



A" arcsin 

2(1 - p) - K - y/K 2 + 4Kp(l - p) 



(18) 



For the downward processes, S(p; K) is replaced by §(1 — 
p;K). 

Out of equilibrium, S is different from the system en- 
tropy S(p) — — plnp — (1 — p)ln(l — p). Indeed, AS is 
the entropy flow and is related to the system entropy 
change AS = S(p(r)) — S(p(0)) via the always-positive 
entropy production AiS — AS — AS > 0. It is only in 
the quasi-static limit, where K and thus AiS = 0, 
that S(p; K) reduces to S(p). 

We are now ready to apply the above results to our 
heat engine. To make the connection with the left/right 
symmetry required for the universality of the coefficient 
in the quadratic term, cf. the discussion in the introduc- 
tion, it will be of interest to consider an asymmetry in 
the rate constant: we will use the subscripts C c and Ch 
for the rate constant C when in contact with the cold and 
hot reservoir, respectively. Recalling that processes I and 
III are upward and downward processes, respectively, the 
boundary condition (ITB)) leads to 

C c t c = F{ Pl -K c ) - F(p ; K c ) (19a) 
C h r h = F(l - p ; K h ) - F(l - pi;K h ), (19b) 

which determine the integration constants K c and Kh, 
respectively. 

Substituting K c and Kh into Eq. (fTT|) . we obtain the 
amount of heat that enters the system during the pro- 
cesses I and III: 

Q c = T c [S( Pl ;K c ) - §( Po ;K c )] - T C AS C (20a) 

Q h = T h [S(l - p ; K h ) - 8(1 - pi;K h ) - T h AS h , 

(20b) 



rj = 1 + 



T C A$ C 
T h AS h ' 



(21) 



In the quasi-static limit, K c — > and Kh — > 0, one has 
S(p;0) = S(p) = S(l-p), hence AS C = -AS h , so that 
Eq. (f2~Tj) reduces to Carnot efficiency. 

The above results provide the required optimization 
with respect to the schedules. It remains to perform the 
optimization with respect to the remaining degrees of 
freedom r c , th, Po and Pi- In general, this can only be 
done numerically since Eq. (|14[) only provides an im- 
plicit equation for the time-dependence of the optimal 
schedule. We are, however, mainly interested in the ver- 
ification of universal features of efficiency at maximum 
power. We therefore proceed with a perturbative analy- 
sis for which analytic solutions can be obtained. 



IV. WEAK DISSIPATION LIMIT 

The deviation from Carnot efficiency can be investi- 
gated using the theory of linear irreversible thermody- 
namics where Th — T c is assumed to be smaller than the 
temperatures Th and T c of the reservoirs. However, for 
finite time thermodynamics, a different kind of expan- 
sion, directly related to the irreversibility caused by fi- 
nite operation time, is more natural. As mentioned in 
the previous section, K is a direct measure of the devi- 
ation from the quasi-static limit. Hence, it is natural to 
expand thermodynamic quantities in K . Since this is an 
expansion about the reversible case of zero dissipation, 
we will refer to this as the limit of weak dissipation. 

We expand Eq. (fl~5|) in a series in \JK. The leading 
term is 



F{p-K) 



arcsin(l — 2p) 



K 



(22) 



With this approximation, we are able to solve Eq. 
for K to obtain 



(a = c, h) 



where 



<pi = arcsin(l — 2pj), (i = 0, 1). 



(23) 



(24) 



The present expansion is thus valid under the following 
condition of weak dissipation: 



C a r a » |0i - </>o|, (a = c,h). 



(25) 



Note that it can easily be satisfied in our model since the 
right hand side is bounded by ir. 

Once we find the value of K, the remaining calculation 
is straightforward. Equation (| 14[) leads to the optimal 
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protocols: 

Pc(t) 

Ph(t) 



1 - sin — |0i - O | - 0o 



1 + sin — |0i - 0o I + 0i 



(26a) 
(26b) 



Expanding in a Taylor series with respect to y K, Eq. 
151) is approximated by the two lowest order terms as 



(27) 



S(p) = 5(p) - arcsin(l - 2p)VK 



Inserting the value of K, we obtain the maximum heat 



= -T C AS-T C 



Ql 



T h AS - T h 



1 ~ 0O) 2 
C C T C 
(01 ~ 0O) 2 
ChTh 



(28a) 
(28b) 



where AS 1 = S(j)q) — S(pi) is the reversible entropy 
change. The second term on the right hand side is the 
irreversible heat, which has to be small under the con- 
dition (|25l) of weak dissipation. In the quasi-static limit 
(t — > oo), the second term vanishes and the efficiency 
([21~j) reaches the Carnot efficiency, as expected. 

When the operation time is too short, the irreversible 
heat becomes dominant and the net heat becomes nega- 
tive. Equation (|28|) indicates that positive power can be 
obtained only if 



(T h - T C )AS T ( 



)) 2 



> 



T h 

ChTh 



(29) 



This inequality is consistent with the condition of the 
asymptotic expansion ()25|) and can thus be satisfied even 
for a large temperature difference. 

So far, we have maximized the power only for the fixed 
operation times r c and Th and the boundary values po 
and pi of the occupation probabilities. Now we further 
maximize the power 



P = 



Qc + Qh 

T c + T h 



(30) 



with respect to the operation times. It is easy to find 
that the power is a maximum when 



M 2 T C (1 + y/T h C c /T c C h ) 
C c AS(T h - T c ) 



0o) 2 ^(l + y/T c C h /T h C c ) 
C h AS(T h - T c ) 



(31a) 



(31b) 



This optimization reflects the usual competition with 
the denominator of the power preferring faster operation 
whereas the numerator suggests a slower schedule to stay 
closer to Carnot efficiency. 

For the asymptotic expansion to be valid, the optimal 
operation times must satisfy the condition (|25p of weak 



dissipation. That is, for the process III the following 
inequality must be satisfied: 



101 - 00 I 



AS 



T h - T 



C, 



hTh 



<h\T h (l + y/T c C h /T h C c ) 



< 1. 



(32) 

This can be achieved in two ways. The first one cor- 
responds to the usual condition for linear irreversible 
thermodynamics, (Th — T c )/Th -C 1. The alternative 
is | AS/ (0i — 0o)| *C 1. In this limit, our result remains 
valid even for large temperature differences. 

With the optimized operation times (|3"Tj). the resulting 
power is written as a function of po and p\ : 



P 



(T h - T c ) 



4(Vn/Ch~ + ^/Tjc~ c y 



: D(p , Pl ), (33) 



where 



D(po,Pi) 



AS 2 



(01 - 0O) 5 



(34) 



The power reaches its maximum when D(po,pi) takes a 
maximum value, D max — 0.439 at Po = Pi = 0.0832 or 
Po = Pi = 0.9168. At these conditions, the optimal oper- 
ation time (|31l) and the maximum heat (I28[) both vanish. 
However, the power remains finite. This final optimiza- 
tion thus leads to a singular and unrealistic situation. 
We note, however, that since Eq. (f34|) does not depend 
on the system parameters, the efficiency does in fact not 
depend on this final optimization step. Therefore, we 
proceed to evaluate efficiency without further reference 
to optimal occupation probabilities. 

Using the maximum heat (|28p and the optimal time 
(f3"Tj) . we finally obtain the following remarkable result 
for the efficiency at the maximum power: 



7? c (i + ychT/CcTh) 

2(1 + y/C h T c /C c T h ) - Vc 



4(1 + y/f) 8(1 + Vr) 



o(rit) (35) 



with r = Ch/C c . When r = 1, the efficiency ([33)) exactly 
coincides with the Curzon-Ahlborn efficiency rjcA = 1 — 
y/T c /Th. Note also that the efficiency is bounded below 
by rj c /2 for C'h/C c — > oo and bounded above by rj c /(2 — 
rjc) for Ch/Cc — > 0. These limits can be realized without 
violating the condition of weak dissipation. 



V. DISCUSSION 

We have calculated the efficiency rj at maximum power 
of a Carnot cycle with a single level quantum dot as 
the operational device. Our calculation is in agreement 
with known universality properties. In particular, the ef- 
ficiency at maximum power is equal to half of the Carnot 
efficiency in the regime of linear response, rj = rj c /2 + .... 
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In the case of a left/right symmetry, corresponding to 
equal exchange rate coefficients Ch = C c of the dot with 
the heat reservoirs, the coefficient of the quadratic term is 
also given by its universal value 1/8, r\ = r/ c /2 + ^/8 + .... 
However, we need to stress that this result was obtained 
not by an expansion in rj c but in the limit of weak dissi- 
pation. In fact, this calculation adds a new perspective 
concerning the occurrence of Curzon-Ahlborn efficiency 
itself. Indeed, in the presence of left/right symmetry, the 
efficiency is actually exactly equal to the CA efficiency 
■qcA = 1 — \/l ~~ Vc m the limit of weak dissipation. This 
limit is reminiscent of the original derivation of CA effi- 
ciency, and is in the present model formally similar to the 
assumption of a linear conduction law between reservoir 
and quantum dot. However the concept of weak dissipa- 



tion is more general. It remains to be explored whether 
this observation implies a wider range of validity of CA 
efficiency. In particular, it could explain why observed ef- 
ficiencies at maximum power are not very different from 
CA efficiency in a wide range of systems under opera- 
tional conditions far from linear response. 
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